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ABSTRACT 

We find that the amount and nature of the assumed ionizing background can strongly affect galaxy 
formation and evolution. Galaxy evolution simulati ons typically incorporate a n ultraviolet background 
which falls off rapidly above z = 3; e.g., that of iHaardt fc Madaul ()1996f ). However, this decline 
may be too steep to fit the WMAP constraints on electron scattering optical depth or observations 
of intermediate redshift (2; ~ 2 — 4) Ly-alpha forest transmission. As an alternative, we present 
simulations of the cosmological formation of individual galaxies with UV backgrounds that decline 
mo re slowly at high redsh i ft: bo th a simple intensity rescaling and the background recently derived 
by iFaucher-Giguere et al.l ()2009D . which softens the spectrum at higher redshifts. We also test an 
approximation of the X-ray background with a similar z-dependence. We find for the test galaxies that 
an increase in either the intensity or hardness of ionizing radiation generically pushes star formation 
towards lower redshifts: although overall star formation in the simulation boxes is reduced by 10 — 25%, 
the galaxies show a factor of ^ 2 increase in the fraction of stars within a 30 kpc radius that are formed 
after z = 1. Other effects include late gas inflows enhanced up to 30 times, stellar half-mass radii 
decreased by up to 30%, central velocity dispersions increased up to 40%, and a strong reduction in 
substructure. The magnitude of the effects depends on the environmental/accretion properties of the 
particular galaxy. 

Subject headings: galaxies: elliptical and lenticular — Galaxy: formation — methods: numerical 



> 

p 
O 

o\ 

O 



X 



1. INTRODUCTION 

The importance of the ionizing background on the gas 
dynamics and hence the star-fo rmation hist o ry of galax- 
ies h as long been recognized ()Rees 1 119861: lEfstathiou l 
I19921) . and studied both n umerically (e.g.lGnedin 1120001) 
and semianalytically (e.g. 'Bens on et al.|[200^ . (In this 
paper we take "ionizing background" to refer to both 
the heating and H/He ionization effects of the 13.6 eV 
- 100 keV portion of the external radiation field expe- 
rienced by galaxies at various epochs; i.e., the energies 
which can ionize H and other species but have a large 
enough cross-section to significantly affect the gas prop- 
erties.) However, although the hydrodynamic simula- 
tion of galaxies has become something of a co ttage in- 
dustry in recent years (see iMaver et al.l l2008l for a re- 
view with special focus on disk galaxies), relatively little 
attention has been paid to the form of ionizing back- 
gr qund which is used. Th e def ault option (seen recently 
in iGovernato et al.l 120071 and iScannapieco et all 120081 
among ot hers) is generally a ver sion of the UV back- 
ground ofl Haardt fc Madai] ()1996[ ) . Those authors mod- 
eled the QSO background and the effects of reprocessing 
by Lya forest clouds. They fit the redshift dependence 
of the H photoionization rate with a generalized Gaus- 
sian of the form (1 + z)^exp [—{z — ZcY / S]\ for their 
parameters {B = 0.73, Zc — 2.3, S = 1.9 in the original 
paper), the intensity declines steeply above z = 3 and 
the background is negligible beyond z = 1 . 

There is, however, a history of simulations using ion- 
izing backgrounds with a less steep redshift dependence. 

Electronic address: |dclayh@astro.princeton.edu| 



iNavarro fc Steinmetz I (|1997l ) is one such example: they 
performed 40'^ SPH particle simulations resamplcd from 
a larger P^M simulation, together with a UV background 
that was constant at high redshift. They found that 
the final amount of cooled gas was reduced by up to 
half, with late-accreted gas preferentially affected, com- 
pared to no background. Moreover, recent evidence such 
as the independent observations of th e optical depth 
to electron scattering (WMAP: Dunkelv et all l2008f ) 
and the Lyman alpha Gun n-Peterson trough (SDSS: 
IChiu. Fan, fc Ostriker |[2003h have put increasingly tight 
constraints o n the reionization history of the universe. 
In particular, iFaucher-Giguere et al.l (|2008( ) examine the 
Lya effective optical depth using 86 quasar spectra, and 
find an essentially fiat ionization rate out to z — 4.2 
(see Fig. [T]and the discussion in §2.1 below); these au- 
thors have thus proposed a new UV background, de- 
scribed in Faucher-Giguere et al. (2009*), which we inves- 
tigate here. DairAgiio et al. (2009), using 1733 quasars 
from SDSS, similarly find a flat ionization rate for 2 < 
z < 4.2, albeit wi t h a sl ightly different normalization. 
(jHaardt fc Madaul ()2001[ ) themselves presented a UV 
background very similar to our New UV below, but that 
model has seemingly failed to gain widespread accep- 
tance.) 

Further, although the diffuse X-ray background is 
becoming increasin gly well constra ined, at least in 
the local universe (iGilh et aD 120071 ) . and its unique 
role in heating the univers e prior to full reion- 
ization has been established ([Venkatesan et al.l l2001l : 
iRicotti Ostriker fc GnedhH l2005f l. few recent simula- 
tions (|Ricotti et al.l l2008l for example) have incorpo- 
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rated this component. The hard X-ray background 
is largely produced by QSOs and lower-luminosity 
AGN. Combining the the total background at z = 
with the observed X-ray spectra of individual sources 
an d the redshift dependence of AGN output permit- 
ted [S^onov^iOitriker^5ISu^ai3 ()2004[ ) to estimate the 
10° — 10^ keV background as a function of redshift. This 
background has a Compton temperature of lO'''^ K due 
to the peak in vJ^ around 30 — 50 keV, which penetrates 
regions that are optically thick to UV (although many 
simulations, including ours, treat the UV as optically 
thin also) and can provide a consid erable source of both 
ioniza tion and heating. Indeed, iMadau fc Efstathiod 
(|1999| ) found that including an X-ray background in- 
creases the equilibrium temperature of the IGM by ~ 
20%. 

Finally, since the formation of massive ellipticals where 
star formation has been effectively quenched since z ~ 
1 — 2 is still not well understood and is the subject of 
ongoing studies (fo r a review see the introduction of 
lHo5k ins et al.ll2008f) . a study to examine the extent to 
which results are sensitive to the assumed ionizing radi- 
ation background is warranted. (Naturally, further work 
could also be done with disk galaxies, Lyman-break pop- 
ulations at higher redshift and others.) 

The paper is organized as follows. In ij2l we de- 
scribe the numerical methods and parameters of our 
simulations, in particular the ionizing radiation back- 
grounds which we apply , comprising a recent version of 
iHaardt k, Madaul ([19961), a new rescaled version of the 
UV background which falls off much more slowly at high 
redshift, this new UV background with an additional X- 
ray component, and the more realistic, recently calcu - 
lated UV background of iFaucher-Giguere et al.l (|2009( ) . 
In [J3]we describe the results obtained from those simula- 
tions, in particular the effects of the backgrounds on the 
gas properties, star formation and stellar dynamics. fJH 
is a discussion of the implications of these results when 
taken collectively, and §5 is conclusion. 

2. SIMULATIONS AND BACKGROUNDS 

2.1. Radiation Backgrounds 

As a baseline, we use the IHaardt fc Madaul (fl99l UV 



background as used in e.g. iNaab et al.l ( 2007 ): we call 
this the "Old UV" model. To create an upper bound 
on high-z UV, we keep the same assumed spectral shape 



but set the intens ity to decline as ( 1 



in physi- 



cal units above the lHaardt fc Madaul peak, which occurs 
at z w 2.4; we call this the "New UV" model. For a 
mo re realistic case, we u s e the background calculated 
by iFaucher-Giguere et al.l ()2009l ). which has a similar 
z-dependence of intensity but whose spectrum softens 
markedly for z > 3 as the quasar contribution dies out; 
we call this "FG UV" . The H ionization rates of these 
th ree m odels arc c ompare d in Fig.[Tl along with the data 
of iFau chcr-Gigucr e et al.l ((2008) mentioned above. It is 
apparent that New UV and FG UV provide a better fit 
to the observations than does Old UV, in particular at 
z > 4. 

We also implement an X-ray background. This 
component uses the spe ctr al s hape given in 
ISazonov. Ostriker. fc SunvaevI ()2004f ). which repre- 
sents an average quasar background considering both 
obscured and unobscured sources, and is strongly 




Fig. 1. — Hydrogen photoionization rate F (s ^) vs. z for the Old 
UV (black dotted line), New UV (blue solid line), New UV+X (red 
dashed line), FG UV (green short dash-dot line) and New UV with 
cutoff (magenta long dash-dot line) models. The data are from 
LEauchcr-Gigucre et al. (2008). Note that our X-ray background 
contributes only slightly to H ionization, though it has significant 
heating and He ionization effects (see text). 

peaked around 30 keV in EFe- Th e intensity normal- 
ization comes from iGilli et al.l (|2007f ). who modeled the 
AGN/QSO X-ray background using both deep pencil- 
beam pointings and shallow surveys. This is converted 
into heating and i onization rate s using Cloudy (v07.02, 
last described in i Ferland et al.l I1998I )'. which includes 
photo and Compton heating as well as secondary ioniza- 
tions; the heating rates are then increased by a factor 
of 1.5 to bett e r agr ee with the more recent model of 
ISazonov et al.l ()2005D . Heating rates due to the UV and 
X-ray backgrounds are roughly equal for virialized gas 
at z = 0. The redshift dependence of this background 
is taken to be similar to New UV: intensity in physical 
units increases as (1 -I- z)'^ to z = 2, and declines as 
(1 + z)^^ thereafter (which is admittedly unrealistic; see 
below). We call this the "New UV-hX" model. Notice 
in Fig. [1] that New UV-I-X has only a negligibly higher 
ionization rate than New UV: photons at keV energies 
and above deposit > 99% of their energy as heat in 
a highly ionized medium through e lectron-electron 
colhsions (IShull fc van Steenbergl Il985f ). The X-ray 
background also contributes at the ~ 30% level to the 
higher-energy Hell ionizations. 

As an additional motivation for the revised ionizing 
background, we calculate a simple (homogeneous) semi- 
analytic model of reionization using Cloudy. We cre- 
ate bins at successive epochs and apply the correspond- 
ing background (based on our z-dependence formulas) 
to gas of the corresponding (physical) density (i.e., the 
mean present density scaled as (1 -I- z)"^); Cloudy out- 
puts the electron-scattering and Gunn-Peterson optical 
depths (tos and tqp, respectively). For t^s this method 
is replaced above redshift 10 by an analytic integral of 
the electron density (i.e., Tcs = ct J Uedl over the ap- 
propriate cosmology, with rig given by the balance of 
photoionization with coUisional recombination); FG UV 
is specified to have zero intensity above z = 10, so in 
that case this integral is also. Note that removing 
the assumption of homogeneity and including the effects 
of dumpiness would in general decrease the effective re- 
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combination time (jMadau. Haardt. fc Reed [l999f ). thus 
ieading to a less-ionized universe and decreasing both 
Tqs and the epoch of reionization; this strengthens our 
conclusions below. Moreover this calculation is optically 
thin and in equilibrium, neither of which is the case for 
real reionization; however our Gadget code also assumes 
optical thinness and equilibrium (see below); so these re- 
sults at least show what our simulations would produce 
for Tes and tgp and may be compared to each other, even 
though they are physically oversimplified and not directly 
comparable to observations or more detailed calculations. 

The results are summ arized in Tab l e [1] and 
compared with WMA P -5 (iDunkelv et al.l l2008l ) and 
iChiu. Fan, fc Ostriker I ([20031) (CFO). The Old UV 
model is significantly {2.7a) low compared to WMAP in 
its prediction of t^s, while New UV and New UV-I-X are 
extremely high (6.7 and 9.2a, respectively). The evolu- 
tion of the Gunn-Peterson (Lya) optical depth tgp shows 
a qualitatively similar pattern, with the Old UV and New 
UV models comfortably bracketing the CFO results. FG 
UV agrees very well with both the WMAP Tcs and the 
CFO Tgp values. 

The values of t^s for the New models may seem at first 
too high to give physically relevant results in a simula- 
tion; however, these models exceed the WMAP value 
only due to contributions from z > 10, where it has 
been shown (jPiikstra et al.l [2004) that the ionizing back- 
ground has comparatively little effect. (Again, FG UV 
is specified to reach zero intensity at exactly z = 10, and 
matches the WMAP t^s value very well.) Even though 
the gas may be highly ionized, its cooling time is still 
very short, so the e ntropy and hence the dyn amics are 
not much affected (|Thoul fc Weinberg I \l9M ). There- 
fore we expect our results not to differ so dramatically 
from a WMAP-consistent model as the Tos results would 
suggest. To test this hypothesis, we implement a high-z 
cutoff for the UV of the form [exp(- ( (2:-2 4) /4.7)2), z > 
2.4], analogous to lHaardt fc Madaul (|1996D with param- 
eters chosen to exactly match the WMAP-5 value of t^s, 
and run a test simulation of Galaxy A to z = 0.5 (this 
model is called "New UV with cutoff" in Tabled. We 
find this cutoff causes essentially no change in gas prop- 
erties for all z < 5.2, and only a small (1.5%) increase 
in total stellar mass at z = 0.5 (by contrast the Old UV 
simulation has 5% more stellar mass than New UV at this 
redshift). As we wish to ensure both that the new models 
represent a firm upper bound, and that any differences 
from Old UV are large enough to be easily discernible, we 
use the models without a cutoff in the rest of our analy- 
sis. These cutoff results also indicate that while New UV 
and FG UV have slightly different redshift dependences 
for their intensities, those differences will be negligible 
for our purposes compared to the differences in spectral 
shape of the two models. In other words, comparing the 
FG UV and New UV results will isolate the effect of the 
spectral shape of the UV background, while comparing 
New UV to Old UV will isolate the effect of its intensity 
(represented by e.g. H photoionization rate) as a func- 
tion of z. In fact, our results will show that for most of 
the properties we study, increasing the intensity and in- 
creasing the spectral hardness give the same qualitative 
effects. 

Furthermore, with regard to the New UV-I-X model, 
we don't generically expect UV and X-ray backgrounds 



TABLE 1 

Reionization with various backgrounds 



Name 



z{tgp = 1) z(rGP 



WMAP-5 
CFO 



0.087 ± 0.017 
0.12 ±0.03 



Old UV 0.040 ± 0.002 

New UV w/ cutoflf 0.087 ± 0.011 

FG UV 0.082 ±0.011 

New UV 0.202 ±0.011 

New UV+X 0.243 ± 0.023 



4.21 ± 0.02 5.95 ±0.02 



4.13 ±0.01 
4.09 ±0.02 
4.21 ±0.02 
4.53 ±0.02 
4.77 ±0.02 



4.78 ±0.01 
5.46 ±0.02 
5.90 ±0.02 
6.28 ±0.02 
6.49 ± 0.02 



Note. — "CFO" results a re from IChiu. Fan, fc Ostriker I ||200^ '): 
"WMAP-5" results are from IDunke v et al.l II2008I ) These observa- 
tional results are italicized. 

to have the same z-dependence; in fact, since the X-rays 
primarily come from quasars only rather than quasars 
and stars, and see a much smaller optical depth, a more 
sudden dropoff in the X-ray background intensity above 
the quasar luminosity function peak at z ~ 2 might be 
expected. However, again we are looking for an up- 
per bound and for any effects to be exaggerated for 
easy detection, and our overestimate is relatively mod- 
est: at z = 3.9, our New UV+X model gives a mean- 
density IGM temperature Tq of - 3.3 x 10^ K, 50% 
hig her than the observatio nally-fit value of 2.2 x 10* 
K (jZaldarriaga et al.ll2001h . 

In sum, the FG UV model of the ionizing radiation 
background seems the best of those considered in match- 
ing the observational constraints of H ionization rate, 
electron-scattering optical depth, and Gunn-Peterson op- 
tical depth. In the following sections we will explore the 
effect of this background on galaxy evolution, in compar- 
ison with other backgrounds that are relatively too weak 
(Old UV) and too strong (New UV, New UV+X). 

2.2. Simulations 

We apply the diff erent radiation field s to the SPH simu- 
lations described in lNaab et al.l (|2007f ): see that paper for 
full details. The code is based on GADGET-2, and the 
galaxies are ellipticals selected from a (50 Mpc/h)^ box 
with cosmological (ACDM) initial conditions, and resim- 
ulated in a (10 Mpc/ft,)^ box centered on each galaxy, 
using high-resolution DM and gas/star particles for the 
volume containing all particles which are within 500 kpc 
of the cental galaxy at z = (this is roughly 1.5 Mpc/h 
comoving). Importantly, the simulation does not in- 
clude optical depth effects, in particular the self-shielding 
of dense star-forming regions from the ionizing back- 
grounds, although those regions would be optically thin 
to X-rays regardless. Neither does it include any feed- 
back from supernovae or AGN; although those processes 
are certainly important in many if not most galaxies, 
our object is show the differential effect of changes in 
the ionizing background on the evolution of galaxies, in- 
dependent of feedback effects. (However, it is interest - 
ing to note that recent work by I Pawlik fc Schavj ()2008D 
suggests that feedback and background effects may not 
be independent, and may in fact amplify one another.) 
Star formation is performed at a fixed density threshold 
(Pcrit = 7 X 10~^^ g cm~^, or nn cnt — 0.03 cm~^, as in 
INaab et al.l[200l . 

For easy comparison, we use the same set of initial con- 
di tions that were de signated galaxies/halos A, C, and E 
in lNaab et al.l (|2007( ). and are so designated here as well. 
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All simulations were performed with 100'^ SPH particles 
(corresponding to a gravitational softening length of 0.25 
kpc for the gas and star particles, and twice that for the 
dark matter particles; gas and star particles have masses 
of the order 10® Af©). As will be discussed below, galaxies 
A and C gave the same qualitative results, while galaxy E 
was somewhat different, due to its different merger and 
accretion history. Sin c e Gal axy A was the most well- 
studied in iNaab et al.l ()2007l ). we choose here to focus 
on it, bringing in the other two galaxies where relevant. 
(FG UV, a later addition to the study, was run on galaxy 
A only.) 

Throughout the paper, all distances are physi- 
cal except where noted; the assumed cosmology is 
(17M,^A,^ b /»M, q-s,fe) = (0.3,0.7,0.2,0.86,0.65) as in 
INaab et all pOOTh . 

3. RESULTS 
3.1. Gaseous effects 

We naturally expect an increase in the ionizing back- 
ground to lead to more efficient gas heating at high red- 
shifts, but the consequences on galaxy formation cannot 
easily be predicted. 

Figure [2] shows the temperature distribution of the gas 
in the central 2 Mpc comoving of the galaxy A simula- 
tion (essentially the high resolution region of validity; the 
simulations, being resampled from a larger box, are non- 
periodic) with the three backgrounds, at z = 5.2 (chosen 
to be where the Old UV background is nonzero but still 
substantially lower than New UV / FG UV) and z = 0. 
At z = 5.2 the models are well separated; in fact the 
gas in the Old UV simulation has just finished heating 
from a very cold (lOOK) state (to which it had cooled by 
adiabatic expansion from its initial state ai z — 24; the 
cooling function implemented in the code is primordial 
(H-He only) and cuts off at 10*K), and has just finished 
H/Hel reionization. As expected, the New UV model, 
which has been entirely reionized (including Hell) since 
z 2± 8, has significantly more hot gas, and including X- 
ray heating pushes most of the gas up to temperatures 
significantly above lO'^K. The FG UV model lies in be- 
tween Old UV and New UV, which is expected since it 
has a higher intensity than the former but a softer spec- 
trum than the latter. FG UV has been H/Hel reionized 
since z ~ 8 but has not yet reionized Hell. 

At z = 0, however, the Old UV, New UV, and FG UV 
models are essentially identical; increasing the intensity 
or changing the spectrum of the ionizing background at 
high redshift has no effect at the present (in other words, 
the cooling time for all gas is shorter than 11 Gyr, since it 
has forgotten the extra heat from z > 2 .4 ). T his is consis- 
tent with Mesi nger. Brvan. fc Haimani (j2006f ). who found 
that the effects of a UV background which is suddenly 
turned off begin to dissipate after 0.3 of a Hubble time; 
since we are essentially turning off the extra radiation at 
z = 2.4, we would expect a similar convergence. More- 
over, iHui fc Haimaiil ()2OO30 found that well (> 2 Gyr) 
after reionization, the IGM equilibrium temperature ap- 
proaches a value that depends only on the spectral shape 
and not the intensity of the ionizing background, so even 
if Old UV and New UV had different intensities all the 
way to 2 = we wouldn't expect a significant difference. 

Adding X-rays (to all epochs), on the other hand, sig- 
nificantly heats the gas, especially the coldest gas: the 



mean gas temperature rises 15% from 1.3 x lO'^ K to 
1.5 X 10^ K. The addition of X-rays also produces a 16% 
larger total mass of gas due to reduced star formation in 
small systems, as we will see in the next section. This 
larger reservoir of cool (10^ — 10^ K) gas (cf. the lower 
panel of Fig. [2|) helps to prolong and enhance the epoch 
of star formation in massive galaxies. The warm-hot gas 
mass (WHIM; lO*-^ < T < lO'^ K) is 30% larger in the 
X-ray case than in the cases without X-rays. 

Of course, we are especially interested in the gas which 
has collapsed and virialized in dense halos. Figure [3] 
shows the temperature spectrum of high-density gas 
(p > 200p, where p is the mean baryonic density of the 
universe) at z = for the four backgrounds and the three 
galaxies. (We note again that these simulations include 
no optical depth effects and therefore overestimate the 
UV flux that virialized regions see.) We see that adding 
early UV does not affect the amount of cold (lO^K) dense 
gas, and has uncertain effect on the hot (lO^K) dense 
gas, increasing it slightly in galaxy C, making negligi- 
ble change in galaxy E and decreasing it somewhat in 
galaxy A (although not for FG UV). Adding an X-ray 
background substantially increases the hot dense gas for 
galaxies A and C, while having negligible effect on galaxy 
E. 

We explain these differences by reference t o the merger 
histor ies of the three galaxies. As explored in lNaab et al.l 
()2007( ). galaxy A has a merger of mass ratio 6.5 : 1 at 
z ~ 0.6 (6 Gyr ago), galaxy C has a merger of mass 
ratio 3.5 : 1 at z « 0.8 (7 Gyr ago), while galaxy E has 
no significant merger events after an equal-mass merger 
at z « 1.5 (10 Gyr ago). Since the gas-to-star ratio 
is a strongly declining function of halo mass, especially 
when there is significant ionizing radiation to keep the 
low-density gas in small halos hot, we expect in New 
UV+X for the accretion of smaller halos at later times 
to add more hot gas compared to accreting larger halos 
at earlier times, and the earlier the gas is added to the 
dense central galaxy, the more of it can cool and form 
stars in situ. On the other hand, the mergers in the New 
UV (no X) case involve large amounts of colder gas (see 
next paragraph) , which mingles with the existing gas and 
forms stars quickly, thus paradoxically resulting in less 
hot gas at the present for Halo A. 

Figure 2] shows the accretion rate of gas onto the cen- 
tral physical lOkpc of galaxy A; galaxies C and E had 
qualitatively similar results. We see immediately that 
New UV and New UV-I-X have significantly higher ac- 
cretion compared to Old UV, especially in the last 4Gyr 
(z < 0.3) (and in the no X-ray case, the peak at ~ 8 
Gyr ago, corresponding to the merger event mentioned 
above); FG UV has slightly enhanced accretion. This 
is easily understood as gas that was kept hot at early 
times by the harde r early background o f New UV and 
New UV-hX (as per lHui fc Haimani 12003?) finally cooling 
and flowing inward; in the model with lower UV, much 
of this gas would have formed small stellar systems and 
been unavailable for late inflows. At z = 0, the inflowing 
gas in the New UV model has a mean temperature of 
4 X 10^ K, while for New UV-I-X the mean temperature 
is 3 X 10® K. 

3.2. Stellar effects 
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Fig. 2. — Temperature distribution of gas at z = 5.2 (top) and 
z = (bottom) for simulation A (r < 2Mpc comoving); C and E 
have similar results and are not shown. At z = 5.2 the Old UV, 
New UV and New UV+X models are clearly separated, while at 
z = the addition of UV at early times has essentially no residual 
effect but the X-rays (which are present at all times) significantly 
heat the gas. 

Tlie effects of the different background models on star 
formation fits well with what we saw in the gas. Figure [5] 
shows the star-formation history of the central 2 Mpc of 
simulation A (that is, the formation history of stars in 
the central 2 Mpc at z = 0). As we would expect, the 
early star star formation peak at z ~ 4 is increasingly 
suppressed as we increase the intensity and hardness of 
the ionizing background, from Old UV to FG UV to New 
UV to New UV-I-X (and this carries forth to the total 
stellar mass at the present, which is decreased by 1.4% 
for FG UV, 7% for New UV, and 24% for New UV-^X 
compared to Old UV; see Tabled]). However, the in- 
creased radiation above z = 2.4 in the New UV model 
has little effect on global star formation below that red- 
shift. In contrast. New UV+X suppresses star formation 
until the last 2 Gyr, where there is a modest bump from 
the inflow of late-cooling gas. 

We also study the star-formation effects in a single, 
dense region: viz., the central 30 kpc of galaxy A, pre- 
sented in Fig. [6l Here there is little effect on the initial 
peak, which we expect since this region should have a 
short cooling time due to its overdensity. However, star 
formation below z = 1 is enhanced: 19% of the stellar 
mass within 30kpc is formed after z = 1 in the New UV 
case, compared to 11% for Old UV (and 10% for FG UV). 
The bump in New UV centered at ^7 Gyr ago matches 
well with the merger-event peak in gas accretion we saw 
centered at ~8 Gyr ago in FigH) This results in the New 
UV version of galaxy A actually having 5% more stars 
at the present within the 30 kpc radius compared to Old 
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Fig. 3. — The same as Fig. [2] (r < 2Mpc; z = 0), but for gas 
with p > 200p, for galaxies A (top), C (middle), and E (bottom). 
Adding more UV at early times has little effect, while adding the 
X-rays removes the cold-dense gas while substantially increasing 
the hot dense gas in two of three cases. 



TABLE 2 

Stellar mass results 



IC Name 


Background 


A/H(5kpc) 


Mfl(30kpc) 


Mr(2Mpc) 


A 


Old UV 


9.28 


14.00 


37.45 


A 


FG UV 


f0.09 


14.14 


36.92 


A 


New UV 


fO.23 


14.72 


34.95 


A 


New UV+X 


9.9f 


13.91 


28.52 


C 


Old UV 


8.55 


14.12 


36.34 


C 


New UV 


9.28 


16.18 


35.58 


C 


New UV+X 


f0.81 


16.69 


29.58 


E 


Old UV 


8.24 


12.43 


27.98 


E 


New UV 


8.78 


16.71 


27.75 


E 


New UV+X 


f2.66 


13.78 


25.29 



Note. — Results are the mass of stars within the specified radius 
of the principal halo at z = 0. Masses are W^^Mq. 
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Fig. 4. — Accretion rates of gas onto the central lOkpc (physical), 
for galaxy A. The New UV and New UV+X models have increased 
accretion at late times (LBT< 4Gyr). (Missing data points indi- 
cate zero accretion or net outflow.) 



Fig. 6. — The same as Fig. \5\ but only considering star parti- 
cles which are within 30 kpc of the central halo at z = 0. There 
is increased star formation at low redshift (z < 1) as radiation 
increases: the mean (mass- weighted) age of these stars decreases 
from 11.4Gyr (Old UV) to 10.9Gyr (New UV-fX). 
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Fig. 5. — Star formation history of the central 2 Mpc of the 
Halo A simulation, for the four backgrounds. The high-z peak is 
suppressed as radiation increases (meaning as we move from Old 
UV to FG UV to New UV to New UV+X, increasing the intensity 
and hardness of the ionizing background). 

UV; FG UV has only a 1% increase (again see Table [5]). 
New UV-I-X has a very late burst of star formation 0.5 
Gyr ago which is due to a the infall of warm gas that 
cools in situ, but this is not enough to make up for the 
earlier deficit, and it has slightly less stellar mass within 
30 kpc than Old UV; 15% of the stellar mass is formed 
after z — 1. 

Finally, we convert these star-formation histories 
into present luminositie s via the GALAXEV code of 
iBruzual fc CharlotI ((200l . assuming a Salpeter IMF and 
solar metallicity; bolometric magnitude and color results 
for Galaxy A are presented in Table [31 No dust extinc- 
tion effects were included. In brief, we see that compared 
to Old UV, New UV has only a slight increase in bolo- 
metric luminosity (0.13 mag) and blueness (0.09 mag), 
and FG UV has even less change, while New UV-I-X is 
1.3 mag or 3.25 times brighter, and 0.28 mag bluer in 
U-B, thanks to its late burst of star formation. 

These results, like those for the gas, are dependent 
on the environment of the galaxy. For example, with 



TABLE 3 

Observational Characteristics, Galaxy A 



Model Name 


Bol. Mag. 


B Mag. 


R Mag. 


U-B 


B-V 


Old UV 


-21.91 


-20.16 


-21.72 


0.59 


0.93 


FG UV 


-22.00 


-20.30 


-21.82 


0.55 


0.91 


New UV 


-22.04 


-20.39 


-21.85 


0.50 


0.87 


New UV+X 


-23.19 


-21.96 


-23.14 


0.31 


0.69 



Note. — Results are for a radius of 30 kpc at 2 = 0. 



the Galaxy E simulations the New UV+X model has 
an 10% decrease in total stellar mass compared to Old 
UV, but has 34% more stellar mass in the central 30kpc 
from enhanced gas accretion. (New UV with no X-rays 
produces 1% less total stellar mass and an 11% increase 
within 30 kpc.) There is also a correspondingly stronger 
effect on the mean stellar age for galaxy E: the fraction 
of stellar mass within 30kpc formed after z = 1 increases 
from 6% for Old UV to 35% with New UV+X. See glfor 
discussion of these differences. 

3.3. Dynamical effects 

In addition to affecting the mass and hydrodynamic 
properties of the stars and gas in the simulations, the 
radiation background affects their arrangement; i.e., the 
dynamical properties of the galaxy. The first column of 
Table [2] shows a significant increase in the stellar mass 
within a 5 kpc radius as one moves from the Old UV to 
the New UV to the New UV+X case for all three galaxies 
(and FG UV again lies between Old UV and New UV): 
the mass increase from Old UV to New UV is an average 
of '^10%, and from New UV to New UV+X averages 
^^20%, although this includes a slight decrease for Galaxy 
A. This increase is natural considering the enhanced gas 
flows we saw in Fig. SI once that gas reaches the center 
it must cool and form stars there. Figure [7] shows the 
circular speed — GM{r)/r for the innermost 2 kpc of 
galaxy A at z = 0. We see that the stars become more 
centrally concentrated as one moves up in background 
intensity or hardness (from Old UV to FG UV to New 
UV to New UV+X); this is a general feature for z < 2.5. 
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Fig. 7. — Circular speed for stars (dashed line), gas (dotted 
line), and total (solid line; also includes dark matter), for Halo A 
with the four backgrounds, 2 = 0. Colors for the four background 
models are as before. The peak value and inner slope increase with 
increasing radiation. 

Residual gas also increases in the same way. 

Figure [5] shows the number of discrete (i.e., non- 
substructure) stellar systems with masses > 3 x 
IO^A/q identified from the star particles with the 
HOP halo-finding algorithm (jEisenstadt fc Hut I [1991) 
for the galaxy A simulation (within 2 Mpc comoving' ) 
over time. M any pas t simulatio ns (Efstathiou 199^ 
Quinn. Katz. fc Efstat hloin Il996l : [Thoul fc Weinbcr£] 
1996( ) have shown that photoionization can prevent the 
formation of smaller stellar systems, and indeed we find 
that the number of total stellar systems decreases from 
Old UV to New UV to New UV-^X, although we find 
the effect of increased high-z UV on the number of inde- 
pendent halos doesn't persist past z = 0.5, where merger 
events consume several small halos that have formed in 
the Old UV case. Interestingly, FG UV shows fewer small 
halos than New UV: we see the suppression of halo for- 
mation (1.5 > z > 0.75) combined with late mergers 
{z < 0.5). In all cases star formation in smaller halos is 
preferentially suppressed. A future paper (Hambrick & 
Ostriker 2009, in prep.) will deal in more detail with the 
effect of radiation backgrounds on the mass spectrum of 
galaxies. 

Figure [9] shows the half-mass radius for stars within 
30kpc for the three radiation models with galaxy A. The 
major effect is secular growth from z = 8 to the present, 
combined with the stochasti c effects of major mergers. 
iNaab. Johansson fc Ostrikeil (j2009) found in a detailed 
examination of galaxy A (with the Old UV background, 
and including SN feedback) that minor mergers are pri- 
marily responsible for the factor of > 3 increase in half- 
mass radius from z = 3 to 0. The radiation background 
does have an effect, however: New UV and New UV-I-X 
show substantially (> 25%) smaller radii compared to 
Old UV at all z < 2.5, modulo the intermittent merger 
effects; FG UV has negligible difference from Old UV, 
and in fact has a slightly larger radius at z = 0. 

For z > 3 the picture looks somewhat different: specif- 
ically, in our snapshot at z = 2.9, the New UV model 
for galaxy A has a 13%, and FG UV a 27% larger 
stellar half-mass radius than Old UV, (although New 
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Fig. 8. — Number of discrete stellar systems with masses 
> 3 X IO^Mq identified by HOP for the galaxy A simulation (within 
2 Mpc comoving) over time. New UV suppresses small halos some- 
what for 0.1 < z < 1.5, while New UV+X does so strongly for 
all 2 < 1.5; FG UV shows both the suppression of halo formation 
(1.5 > 2 > 0.75) we see in New UV and the late mergers (2 < 0.5) 
of Old UV, resulting in fewer systems than either (though still 
many more than New UV-I-X. 
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Fig. 9. — Half-mass radii of stars within 30 kpc for Halo A. The 
increased radiation models show consistently smaller radii. (The 
large jumps are the result of merger events.) 

UV-I-X is smaller than Old UV by 13%). Further, 
New UV and New UV-I-X both show a 5% reduction 
in pea k circular speed compared t o Old UV at that red- 
shift. Joung, Gen & Brvan' (200§) found, using a mod- 
ified Haardt & Madau 1996 background similar to our 
"New UV with cutoff" , that simulated galaxies have half- 
light radii that are too small and peak circular speeds 
that are too large compared to observed galaxies at red- 
shift 3. However, these authors also used a more strin- 
gent star-formation criterion: when they repeated their 
simulations using the same density threshold as in this 
work, their results at z = 3 agreed with obesrvations 
(Ryan Joung, priv. comm.). 

We speculate that the increase in central concentration 
with increasing radiation is enhanced by reduced sub- 
structure in the galaxy (because small stellar halos are 
suppressed) and correspondingly less dynamical friction. 
To test this hypothesis, we have constructed a statistic 
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TABLE 4 

Substructure Measure 



Model Name 


A 


C 


E 


Old UV 


114.6 


98.9 


113.9 


FG UV 


19.1 






New UV 


14.2 


13.7 


129.8 


New UV+X 


31.5 


17.8 


70.6 



Note. — Data are cruM(4kpc); see 
text for definition. 



to measure the substructure or "lumpiness" that we des- 
ignate cruM(''), which is created by an unsharp-mask-Hke 
procedure: 

1 ^ 

• Rc/5<R(i)<2He 
U{i)>0 

^^'^ {poGirm ' 

where Re is the effective (half-mass) radius of the galaxy, 
p is the stellar density field (created by cloud-in-cell map- 
ping of the star particles to a grid of resolution 0.38 kpc, 
or twice the gas softening length), and G(r) is a gaus- 
sian of width r. That is, we create a mask by smoothing 
the stellar density field of the galaxy with a fixed-width 
gaussian kernel, then sum ((data — mask) /mask) ^ over 
all pixels in the grid where it is positive (i.e. overdense 
regions) except the central peak, and normalize to the 
size of the galaxy. We find that the Old UV model in- 
deed has far more substructure than the two new models, 
at least for halos A and C; Table 2] shows cruM(4kpc) for 
th e three backgroun ds and ICs. Halo E, as remarked on 
in lNaab et al.l (|2007t l. is composed mainly of stars which 
formed in-situ rather than the accretion of smaller stellar 
systems, and therefore we expect its substructure to be 
much less affected by the ionizing background, which is 
what we find. 

At first glance. Fig. [5] and Table H] may seem to be 
incompatible, since the two models without X-rays are 
nearly the same in the former and vastly different in the 
latter. However, if we increase the outer radius in the def- 
inition of (TuM from two to three effective radii (roughly 
from 20 to 30 kpc for Galaxy A), we find that (TuM(4kpc) 
becomes (125.7,145.8,31.5) for Old UV, New UV, and 
New UV+X respectively, a result much more in accor- 
dance with Fig. [H That is, the New UV simulation has 
numerous small halos, but they are at larger radii from 
the central galaxy than in the Old UV case. On the other 
hand. Old UV and New UV have the same number of in- 
dependent halos at the present, as seen Fig. \8\ but Old 
UV has significantly more substructure. This means that 
Old UV formed more small halos initially, but they were 
accreted onto the central galaxy by z ~ 0.5, enhancing 
its substructure. 

Figure [TO] shows contours of stellar density for the cen- 
tral 20kpc of galaxy A with the four background models. 
The smaller size of New UV-I-X is readily apparent, as is 
the reduction in the number and mass of subhalos. One 
can also see the sharp central peak in New UV and New 
UV+X which is absent in Old UV. 

4. DISCUSSION 



Taken collectively, our results suggest a picture of the 
effects of an increased early UV radiation background 
along the following lines. At early times (z > 2), the 
more intense/harder radiation and earlier reionization 
make the gas hotter, especially the less-dense gas, which 
cannot cool as effectively (Fig. That is to say, the 
processes of H and He reionization create large injec- 
tions of heat, so even though the gas begins to cool one 
reionization is complete, an earlier reionization means 
the gas spends less time in the cold, neutral state it 
has been in since recombination. At very early times 
(z > 4) the gas has not yet equilibrated, and so the 
more intense but softer FG UV background gives a higher 
temperature than Old UV, but by z = 3.2 the situa- 
tion has reversed, consistent with lHui fc HaimanI ()2003D 
. At late times, this means that the primary galaxy 
has less substructure, since there are fewer small stel- 
lar systems to accrete; it contains more gas in its cen- 
tral regions (Fig. |4]), and hence more late in-situ star 
formation (Fig. [6]), since that gas would otherwise have 
formed stars before falling in; it is still smaller but now 
more tightly bo und due to gaining less energy from dy- 
namical friction (jNaab. Johansson fc Ostrikeill2009f ) and 
more central star f ormation from inf ailing cold streams . 
In the terms of the iBirnboim. Dekel. fc Neistein I ()2007D 
picture, the cold gas inflows persist longer when early 
star formation is suppressed by ionizing radiation. 

We can distinguish the effects of intensity and spectral 
shape by comparing New UV and FG UV, since they 
have roughly the same intensity (H photoionization rate, 
etc.) from 2 < z < 8 but FG UV has a significantly 
softer spectrum due to the rapid falloff of the quasar 
contribution. We find that FG UV in general gives re- 
sults intermediate between Old UV and New UV, for gas 
temperature, rate of gas accretion by the central galaxy, 
star-formation rate, and concentration of stellar mass; in 
other words, an increase in the intensity of early radia- 
tion and an increase in the hardness of the background 
spectrum produce qualitatively the same results in most 
of properties we study here. One exception is in the to- 
tal number of independent stellar systems: FG UV has 
somewhat fewer than both Old and New UV (though still 
far more than New UV+X). 

Continuing the investigation of spectrum, the effects of 
X-ray as opposed to UV radiation are equally interest- 
ing. Since adding the X-rays increases radiation levels at 
all z, not just z > 2.4, the gas at the present is substan- 
tially hotter, and star form ation is pushed to even later 
times, in agreement with [Johansson. Naab fc Ostrikerl 
()2009t ). who found that the accretion of small stellar 
clumps (minor mergers) was sufficient to suppress star 
formation after z 1: in the New UV+X case the 
small stellar clumps are themselves suppressed, so the 
central galaxy continues star formation to the present. 
There is less star formation overall, however, and there- 
fore more dense (hot) gas. Substructure is further re- 
duced and compactness increased (again agreeing with 
iNaab. Johansson fc Ostrikeill2009l ). 

These effects are, however, depen dent on the accre- 
tion history of the galaxy in question. iNaab et al.l ()2007f ) 
described two distinct mechanisms by which galaxies as- 
semble their stellar mass: the accretion of existing stellar 
systems, and the in situ formation of stars from inflowing 
gas. In our sample, the former mechanism is dominant 
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Fig. 10.— Contours of stellar density in the central 20kpc of galaxy A at 2 = for the Old UV, FG UV, New UV, and New UV+X 
models. The scales are kpc; the contours are quarter-decade steps in stellar density, with the inner bold contour corresponding to 3Mq/pc^. 
The reduction in size and substructure (i.e. the number and size of subhalos) for New UV models are both visible, as is the increased 
central peakiness. 

for galaxies A and C, and the latter for galaxy E. The 
increase in ionizing radiation creates a bias toward the 
gas-accretion mechanism by suppressing star formation 
in small halos. Thus in the galaxies where stellar accre- 
tion was important in the low-radiation (Old UV) case, 
A and C, we find that the increase of in situ star forma- 
tion is more or less balanced by a decrease in accreted 
stellar mass, so the total stellar mass in the galaxy at 
the present is not strongly affected, although the mean 
stellar age and hence the luminosity and color are. How- 
ever, in galaxies where stellar accretion is not important 
(like galaxy E), there is no mass loss from that source, so 
extra surrounding gas falls in and leads to a higher total 
galaxy stellar mass. That is, gas which in the Old UV 
case formed stars in small satellite halos gets heated up 
at early times, then at late times {z < 1), flows in and 
forms stars in the large central galaxy. 



5. CONCLUSIONS 



As we develop accurate, detailed modeling for galaxy 
formation based on realistic cosmological initial condi- 
tions, we are learning just how sensitive the results are 
to the physical input parameters. One important param- 
eter is the ionizi ng background radiation field. The com- 
monly accepted iHaardt &: Madaul ()1996D model, while it 
has been of great utility, fails to produce a good match 
to the most relevant observational constraints: first, the 
ionization rate as determined by Gunn-Peterson observa- 
tions of high-redshift quasars, and second, the electron 
scattering optical depth as determined by WMAP. 

We have considered new mo dels for the ionizin g back - 
ground: a simple rescaling of IHaardt fc Madaul ()1996f ) 
both with and without a significant X-ray component, 
intended to represent an u pper bound of what is possible , 
and the recent results of iFaucher-Giguere et al.l (|2009f) . 
representing a reasonable middle ground. We have per- 
formed a state-of-the-art set of cosmological simulations 
to assess the sensitivity of the results to the assumed 
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ionizing background. 

We find that the gas properties at late times are much 
more affected by the X-ray component than by early UV 
enhancement, with the result that there is as much as a 
30% increase in the WHIM component (T > lO^-^K), and 
a four-fold increas in the hot-dense component (T w 10^ 
K; p > 200 p) when X-ray heating is present. However 
even our rescaled background without X-rays (New UV) 
reduces the formation of stars in small systems and allow 
cold flows to persist to later times, markedly increasing 
the amount of late-time {z < 1) star formation in massive 
galaxies. Correspondingly there is less substructure in 
the massive systems due to reduced accretion of smaller 
stellar systems and consequently less gravitational heat- 
ing. The systems do still grow in size due to accretion 
of satellite systems but the effect is less pronounced with 
when X-ray heating. Despite the absence of feedback, 
and at a resolution of 100^ particles (with gas soften- 
ing length of 0.1625/h kpc), all simulations have a star- 
formation rate of less than IMq/ji at z = with the 
exception of the model with X-rays, in which cold flows 
persist to the present. Finally, the mass function for 



small-mass systems is somewhat reduced by extra early 
UV and substantially so with the addition of X-rays; this 
is a topic which will be explored in greater detail in a fu- 
ture paper. 

The background of iFaucher -Gigue re et all (j2009) (FG 
UV), in contrast with New UV, matches both observed 
ionization rate and electron-scattering optical depth, and 
reduces the substructure in massive systems without cre- 
ating galaxies that have too much gas, too many new 
stars, and too small a size at late times. Therefore we 
recommend its adoption over all the models studied here 
for future galaxy simulations. 
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